Testing and tuning symplectic integrators for the hybrid Monte Carlo algorithm in lattice QCD.
نویسندگان
چکیده
We examine a new second-order integrator recently found by Omelyan et al. The integration error of the new integrator measured in the root mean square of the energy difference, 1/2, is about 10 times smaller than that of the standard second-order leapfrog (2LF) integrator. As a result, the step size of the new integrator can be made about three times larger. Taking into account a factor 2 increase in cost, the new integrator is about 50% more efficient than the 2LF integrator. Integrating over positions first, then momenta, is slightly more advantageous than the reverse. Further parameter tuning is possible. We find that the optimal parameter for the new integrator is slightly different from the value obtained by Omelyan et al, and depends on the simulation parameters. This integrator could also be advantageous for the Trotter-Suzuki decomposition in quantum Monte Carlo.
منابع مشابه
Testing and tuning new symplectic integrators for Hybrid Monte Carlo algorithm in lattice QCD
We examine a new 2nd order integrator recently found by Omelyan et al. The integration error of the new integrator measured in the root mean square of the energy difference, 〈∆H2〉1/2, is about 10 times smaller than that of the standard 2nd order leapfrog (2LF) integrator. As a result, the step size of the new integrator can be made about three times larger. Taking into account a factor 2 increa...
متن کاملAlgorithms for Dynamical Fermions
This is the write-up of three lectures on algorithms for dynamical fermions that were given at the ILFTN workshop ‘Perspectives in Lattice QCD’ in Nara during November 2005. The first lecture is on the fundamentals of Markov Chain Monte Carlo methods and introduces the Hybrid Monte Carlo (HMC) algorithm and symplectic integrators; the second lecture covers topics in approximation theory and the...
متن کاملInstabilities in Molecular Dynamics Integrators used in Hybrid Monte Carlo Simulations
We discuss an instability in the leapfrog integration algorithm, widely used in current Hybrid Monte Carlo (HMC) simulations of lattice QCD. We demonstrate the instability in the simple harmonic oscillator (SHO) system where it is manifest. We demonstrate the instability in HMC simulations of lattice QCD with dynamical Wilson-Clover fermions and discuss implications for future simulations of la...
متن کاملHigher Order Hybrid Monte Carlo at Finite Temperature
The standard hybrid Monte Carlo algorithm uses the second order integrator at the molecular dynamics step. This choice of the integrator is not always the best. We study the performance of the hybrid Monte Carlo algorithm for lattice QCD with higher order integrators in both zero and finite temperature phases and find that in the finite temperature phase the performance of the algorithm can be ...
متن کاملChoice of Integrator in the Hybrid Monte Carlo Algorithm
We study efficiency of higher order integrator schemes for the hybrid Monte Carlo (HMC) algorithm. Numerical tests are performed for Quantum Chromo Dynamics (QCD) with two flavors of Wilson fermions. We compare 2nd, 4th and 6th order integrators at various quark masses. The performance depends on both volume and quark mass. On currently accessible large lattices ( V ∼ 24 ), higher order integra...
متن کاملذخیره در منابع من
با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید
برای دانلود متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید
ثبت ناماگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید
ورودعنوان ژورنال:
- Physical review. E, Statistical, nonlinear, and soft matter physics
دوره 73 3 Pt 2 شماره
صفحات -
تاریخ انتشار 2006